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1 Introduction 

There is a strange phenomenon about the generating functions that count lattice walks restricted to the quarter 
plane: depending on the choice of the set 6 C {^/, , \, |, Z', \, |} of admissible steps, the generating 
function is sometimes rational, sometimes algebraic [but not rational], sometimes D-finite [but not algebraic], 
and sometimes not even D-finite. This is quite in contrast to the corresponding problem in ID, where the 
generating functions invariably are algebraic [^. Much progress was made recently on understanding why 
this is so, and only very recently, Bousquet-Melou and Mishna [^ have announced a classification of all the 
256 possible step sets into algebraic, transcendental D-finite, and non-D-finite cases, together with proofs for 
the algebraic and D-finite cases and strong evidence supporting the conjectured non-D-finiteness of the others. 

As usual, a power series S{t) E Q[[t]] is called algebraic if there exists a bivariate polynomial P{T, t) in 
Q[T, t] such that P{S{t), t) — 0, and transcendental otherwise. Also as usual, a power series S{t) is called 
D-finite if it satisfies a linear differential equation with polynomial coefficients. (Every algebraic power series 
is D-finite, but not vice versa.) At first glance, it might seem easy to prove that a power series is algebraic 
or D-finite: just come up with an appropriate equation, and then verify that the series satisfies this equation. 
But as far as lattice walks are concerned, most proofs given so far are indirect in that they avoid exhibiting the 
equation explicitly but merely are satisfied showing its existence. This is probably so because the equations 
appearing in this context are often too big to be dealt with by hand. 

Nevertheless, it is interesting to know the equations explicitly, because they provide a standard canonical 
representation for a series, from which lots of further information can be extracted in a straightforward manner. 
By applying a well-known technique from computer algebra (in modern fashion, cf. Section ||), we have 
systematically searched for differential equations and algebraic equations that the series counting the walks in 
the quarter plane satisfy. These are given in Section]^. We have also made a first step towards classifying walks 
in I? confined to the first octant (cf. Section ^) by considering all step sets 6 with up to five elements, and 
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performed a systematic search for equations of the corresponding series. More than 2000 hours of computation 
time have been spent in order to analyze about 3500 different sequences. 

We do not provide proofs that the equations we found are indeed correct, but the computational evidence in 
favor of our equations is striking. We have no doubt that all the equations we found are correct. In principle, 
it would be possible to supplement the "automatically guessed" equations by computer proofs in a systematic 
fashion, using techniques that have recently been applied to some special cases p5| , p4| , |||. But we found that 
the computational cost for performing these automated proofs would be by far higher than what was needed 
for the mere discovery. 

2 Methodology 

To study generating functions for lattice walks, we follow a classical scheme in experimental mathematics. 
It is based on the following steps: (SI) computation of high order expansions of generating power series; 
(S2) guessing differential and/or algebraic equations satisfied by those power series; (S3) empirical certifi- 
cation of the guessed equations (sieving by inspection of their analytic, algebraic and arithmetic properties); 
(S4) rigorous proof, based on (exact) polynomial computations. 

In what follows, we only explain Steps (SI), (S2) and (S3). A full description of Step (S4) is given in [|[]. 
By way of illustration, we choose an example requiring computations with human-sized outputs, namely the 
classical case, initially considered by Kreweras [ p7[ ^ |||, of walks in the quarter plane restricted to the step 
sete = {^,/,i}. 

2. 1 Basic Definitions and Facts 

We focus on 2D and 3D lattice walks. The 2D walks that we consider are confined to the quarter plane N^, 
they join the origin of to an arbitrary point (i, j) e N'^, and are restricted to a fixed subset S of the step set 
{y, \, t, y, \, I}. If /(n; denotes the number of such walks of length n (i.e., using n steps 
chosen from 6), the sequence /(n; i, j) satisfies the multivariate recurrence with constant coefficients 

/(n + l;i,j)= ^ f{n-i-h,j-k) for n,i,j>0. (1) 

{h,k)ee 

Together with the appropriate boundary conditions 

/(0;0,0) = 1 and f{n;i,j) = if f < or j < orn < 0, 

the recurrence relation ([!]) uniquely determines the sequence /(n; i, j). As is customary in combinatorics, we 
let 

n>0 i,j>0 

be the trivariate generating power series of the sequence f{n;i,j). As f{n;i,j) = as soon as i > n or 
j > n, the inner sum is actually finite, and so we may regard F{t; x, y) as a formal power series in t with 
polynomial coefficients in Q[a;, y]. 

Specializing F{t; x, y) to selected values of x and y leads to various combinatorial interpretations. Set- 
ting X = y = 1 yields the power series F{t\ 1, 1) whose coefficients count the total number of walks with 
prescribed number of steps (and arbitrary endpoint); the choice x = y = Q gives the series F{t\ 0, 0) whose 
coefficients count the number of walks returning to the origin; setting x — 1, y = Q yields the power series 
whose coefficients count the number of walks ending somewhere on the horizontal axis, etc. 

By [|o[ Th. 7], multivariate sequences that satisfy recurrences with constant coefficients have moderate 
growth, and thus their generating series are analytic at the origin. The next theorem refines this result in our 
context. 
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Theorem 1 The following inequality holds 

f{n;i,j)<\&r for all {i,j,n)^n\ (2) 

In particular, the power series F{t; 0, 0), F{t; 1, 0), F{t; 0, 1) and F(t; 1, 1) are convergent in C[[t\] at t — 
and their radius of convergence is at least 1/ 1 6 1. 



Proof: The total number of unrestricted n-step walks starting from the origin is |6|", so the number of walks 
restricted to a certain region is bounded by this quantity. This implies that the coefficient of in F{t; 1, 1) is 
at most |(S|". The bound also applies to the coefficient of in F{t; a, (3) for a, G {0, 1}, as these series 
count walks which are subject to further restrictions. □ 



2. 1. 1 D-finite generating series of walks are G-functions 

A power series Sit) = E„>o i" ^[M] is called a G-functionQ if (a) it is D-finite; (b) its radius of 
convergence in C[[t]] is positive; (c) there exists a constant C > such that for all n e N, the common 
denominator of oq, . . . , a„ is bounded by C". 

Examples of G-functions are the power series expansions at the origin of log(l — t) and (1 — t)" for a G Q. 
More generally, the Gauss hypergeometric series 2^1 (a, 7; t) with rational parameters a, (3, 7, is also a G- 
series [p^. A celebrated theorem of Eisenstein assures that any algebraic power series must be a G-function 
(if S is algebraic, there exists an integer G G N such that anC^'^^ is an integer for all n.) The fact that 
G-functions arise frequently in combinatorics was recently pointed out by Garoufalidis [|o[|. 

G-functions enjoy many remarkable properties. Chudnovsky prove d that the minimal order differential 



equation satisfied by a G-series must be globally nilpotent (see Section |2.4.4| below for the definition and 
an algorithmic use of this notion). By a theorem of Katz and Honda [^2\ [2l|, the global nilpotence of a 
differential operator implies that all of its singular points are regular singular points with rational exponents. 
See also ITsj |l7| for more details on this topic. 

Theorem 2 Let S{t) be one of the power series F{t; 0,0), F{t; 1,0), F{t; 0, 1) and F{t; 1, 1). If S is D-finite, 
then S is a G-series. In particular, its minimal order homogeneous linear differential equation is Fuchsian 
and it has only rational exponents. Moreover, the coefficient sequence of S{t) is asymptotically equivalent to 
a sum of terms of the form Kp^n°'{\ogn)^ for some constants k G K, a G Q, p G Q, and /? G N. 



Proof: The conditions (a) and (c) in the definition of a G-function are clearly satisfied. The only non-trivial 
point is the fact that the series S has a positive radius of convergence in C. This follows from Theorem |l]. 
The Fuchsianity of the minimal equation for S, and the rationality of its exponents, follow by combining the 
results by Katz, Honda and Chudnovsky cited above. The claim on the asymptotics of the coefficients of S{t) 
is a consequence of [ p0[ Prop. 2.5]. □ 

For 3D walks, the definitions are analogous. The trivariate power series F(t; x, y) is simply replaced by 
the generating series G(t; x,y,z) G Q[x,y, z][[t]] of the sequence g{n; i, j, k) that counts walks in N'^ starting 
at (0, 0, 0) and ending at (i, j, k) G N'^. Note that the appropriate versions of Theorems [l] and || hold; in 
particular, the generating series of octant walks G{t; 1, 1, 1) is a G-series whenever it is D-finite. 

2.2 Computing large series expansions 

The recurrence ([l]) can be used to determine the value of f{n; for specific integers n,i, j G N. Theorem ^ 
implies that f{n; is a non-negative integer whose bit size is at most 0{n). If G N, the values f{n; 



The usual definition is more general, the coefficients of S can be taken in an arbitrary algebraic number field. For our purposes it is 
sufficient and convenient to restrict to rational coefficients. 
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for < n,i,i < N can thus be computed altogether by a straightforward algorithm that uses 0{N^) arith- 
metic operations and 0{N'^) bit operations. (We assume that two integers of bit-size N can be multiplied in 
0{N) bit operations; here, the soft-O notation 0{ ) hides logarithmic factors.) The memory storage require- 
ment is proportional to N^. The same is also true for the truncated power series F/v = F{t; x, y) mod . 
For our experiments in 2D, we have chosen N = 1000. With this choice, the computation of the f{n; i,j) is 
the step which consumes by far the most computation time in our calculations.p^ 

Example 1 The Kreweras walks satisfy the recurrence 

f{n + l,i,j) = f{n,i + + f{n,i,j + 1) + f{n,i- l,j - 1) for n,i,j > 0, 

which allows the computation of the first terms of the series F(t; x, y) 

F{t; x,y)^l + xyt + {x^y^ +y + x)t^ + {x^y^ + 2xy^ + 2x^y + 2)t^ 
+ [x^y^ + 3x2y3 + 3a;^y2 + 2y'^ + Qxy + 2x'^)t^ 

+ {x^y^ + 4x^y** + Ax'^y'^ + 5xy^ + 12x^y'^ + 5x^y + 8y + 8x)t^ H 

and also the first terms of the generating series F{t; 1, 1) for the total number of Kreweras walks 

F{t; 1,1) ^ l+t + 3t^ + 7t^ + nf^ + + 12bf + 333t^ + 939t^ + 2597^*^-^ 
7183i^" + 20505^" + 57859*^^ + 163201i" + 469795^^"* + • • • 

In the 3D case, the values g{n;i,j,k) for < n,i,j,k < N can be computed in 0{N'^) arithmetic 
operations, 0{N^) bit operations and 0{N'^) memory space. In practice, we found that computing G mod 
with N = 400 is feasible. 

2.3 Guessing 

Once the first terms of a power series are determined, our approach is to search systematically for candidates 
of linear differential equations or of algebraic equations which the series may possibly satisfy. This technique 
is classical in computer algebra and mathematical physics, see for example [ pi] , |T], ^s]]. Differential and 
algebraic guessing procedures are available in some computer algebra systems like Maple and Mathematica. 

2. 3. 1 Differential guessing 

If the first N terms of a power series S G Q[[t]] are available, one can search for a differential equation 
satisfied by S at precision N, that is, for an element C in the Weyl algebra Q[t]{Dt) of differential operators 
in the derivation Dt — ^i'^^ polynomial coefficients in t, such that 

£(5') = Cr{t)S^''\t) + ■■■ + ci{t)S'{t) + CQ{t)S{t) ^ mod i^. (3) 

Here, the coefficients co{t),.. . ,Cr{t) E Q[t] are not simultaneously zero, and their degrees are bounded 
by a prescribed integer d > 0. By a simple linear algebra argument, if d and r are chosen such that {d + 
l)(r + 1) > iV, then such a differential equation always exists. On the other side, if d, r and N are such that 
{d + l){r + 1) ^ N, the equation translates into a highly over-determined linear system, so it has no 
reason to possess a non-trivial solution. 

The idea is that if the given power series S{t) happens to be D-finite, then for a sufficiently large N, a 
differential equation of type (thus satisfied a priori only at precision N) will provide a differential equation 

We have carried out our computations on various different machines whose main memory ranges from 8 Gb to 32 Gb and which are 
equipped with (multiple) processors all running at about 3GHz. 
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which is really satisfied by S{t) in Q[[t]] (i.e., at precision infinity). In other words, the D-finiteness of a power 
series can be (conjecturally) recognized using a finite amount of information. 

Given the values d, r, N , and the first N terms of the series S, a candidate differential equation of type ^ 
for S can be computed by Gaussian elimination in 0{N^) arithmetic operations and 0{N^) bit operations. 
Actually, a modular approach is preferred to a direct Gaussian elimination over Q. Precisely, the linear algebra 
step is performed modulo several primes p, and the results (differential operators modulo p) are recombined 
over Q via rational reconstruction based on an effective version of the Chinese remainder theorem. (See [ p3[ ] 
for an implementation of this technique in Mathematica.) 

If no differential equation is found, this definitely rules out the possibility that a differential equation of 
order r and degree d exists. This does not, however, imply that the series at hand is not D-finite. It may 
still be that the series satisfies a differential equation of order higher than r or an equation with polynomial 
coefficients of degree exceeding d. 

Asymptotically more efficient guessing algorithms exist, based on fast Hermite-Pade approximation [Q] of 
the vector of (truncated) power series [S, S' , . . . , S*^''-*]; they have arithmetic complexity quadratic or even 
softly-linear in N. Such sophisticated algorithms were not needed to obtain the results of this paper, but they 
have provided crucial help in the treatment of examples of critical sizes (e.g. guessing with higher values of 
d, r, N and/or over a parametric base field like Q{x) instead of Q) needed for the proof in |^]. 

Example 2 (continued) N ~ 100 tenns of the generating series F{t; 1,1) of the total number of Krew- 
eras walks are sufficient to conjecture that F{t; 1, 1) is D-finite, since it verifies the differential equation 
Ci.i{F{t; 1, 1)) = mod , where 

Cis = 4t2(< + l)(3t - 4)(3i - lf{9f + 3t + l)Dj 

+ 2i(3i - l)2(2916i^ - 1296t^ - ib&U^ - mt^ - 93i + h2)Dl 

+ 3(3i - l)(29808t^ - 26244t^ - 28440i'' + 27Mt^ + 43^^ + 448t - 'm)Df (4) 
+ 6(68040i^ - 88452i^ - 37206t^ + 16758t^ + 954t2 + 253i - 126)A 
+ 18(6480i^ - 8856i^ - 3078t^ + ll'if + 211< + 2). 
Thus, with high probability, F{t; 1, 1) verifies the differential equation Ci,i{F{t] 1, 1)) — 0. 

Sometimes (see Section 2.4.4) one needs to guess the minimal-order differential equation C„iin{S) — 
satisfied by the given generating power series. Most of the time, the choice (c?, r) of the target degree and 
order does not lead to this minimal operator Worse, it may even happen that the number of initial terms N 
is not large enough to allow the recovery of £min, while these N terms suffice to guess non-minimal order 
operators. (The explanation of why such a situation occurs systematically was given in for the case of 
differential equations satisfied by algebraic functions.) A good heuristic is to compute several non-minimal 
operators and to take their greatest common right divisor; generically, the result is exactly £inin- 

As a final general remark, let us point out that a power series satisfies a linear differential equation if and only 
if its coefficients satisfy a linear recurrence equation with polynomial coefficients. A recurrence equation can 
be computed either from a differential equation, or it can be guessed from scratch by proceeding analogously 
as described above for differential equations. 

Example 3 (continued) N = 100 terms of the series S{t) = F{t] 1, 1) suffice to guess that its coefficients 
satisfy the order-6 recurrence 

2(n + 6)(n + 7)(2n + 13) (7n + 34)u„+6 - {n + 6)(140n^ + 2402n^ + 13687n + 25843)u„+5 
+ 3(28n'* + 62671^ + 512371^ + 18281?! + 24070)u„+4 

- 18(n + 4)(28n3 + 311n^ + 897n + 304)m„+3 + 108(n + 3)(35n^ + U3n^ + 1787n + 2309)u„+2 

- 324(n + 2)(7n^ + 90n^ + 382n + 545)u„+i - 972(n + l)(n + 2){n + 4)(7n + 41)u„ = 0. 
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2.3.2 Algebraic guessing 

If the first N terms of a power series S G Q[[t]] are available, one can also search for an algebraic equation 
satisfied by S at precision N, that is, for a bivariate polynomial P{T, t) in Q[T, t] such that 

P{S{t),t)^Cr{t)S{ty + --- + ci{t)S{t)+co{t)^Qmodt^. (5) 

A similar discussion shows that candidate algebraic equations of type for S can be "guessed" by performing 
either Gaussian elimination or Hermite-Pade approximation on the vector [1, S*, . . . , S^], followed by a gcd 
computation in Q[T, t] applied to two (or more) different guesses. 

Example 4 (continued) N — 100 terms of the series S{t) = F{t; 1, 1) counting the total number of Kreweras 
walks suffice to guess that F(t; 1, 1) is very probably algebraic, namely solution of the bivariate polynomial 

Pui{T,t) = t'iit - ifT^ + 6t^(3t - l)^r^ + t^iit - l)(135t^ - m + U)T'^ 

+ 4i2(3t - l)(45i2 - m + 4)r3 + t(3t - - 2Qt + 9)T^ (6) 

+ 2{3t ~ 1)(27<2 - 2i + 1)T + 43t2 + i + 2. 

2. 4 Empirical certification of guesses 

Once discovered a differential equation (||) or an algebraic equation (||) that the power series S{t) seems 
to satisfy, we inspect several properties of these equations, in order to provide more convincing evidence 
that they are correct. These properties have various natures: some are computational features (moderate bit 
sizes), others are algebraic, analytic and even arithmetic properties. We check them systematically on all the 
candidates; if they are verified, as in the Kreweras example, this offers striking evidence that the guessed 
equations are not artefacts. 

2.4. 1 Size sieve: Reasonable bit size 

The differential equation (^ has typically much lower bit size than a differential equation produced by the 
same guessing procedure applied to the same order, degree and precision, but to an arbitrary series having 
coefficients of bit-size comparable to that of S(t). A similar observation holds for the algebraic equation (^. 

Example 5 (continued) If we perturb the coefficients of S{t) = F(t; 1, 1) by just adding a random integer 
between —100 and 100 to each of its coefficients, then the differential guessing procedures at order r — A, 
degree d = 9 and precision N = 100 will either give no result (the over-determined system approach) or 
produce fake candidates (the Hermite-Pade approach) with polynomial coefficients in t, whose coefficients 
in Q have numerators and denominators of about 500 decimal digits each, instead of 4 digits for Ci^i. 

2.4.2 Algebraic sieve: High order series matching 

The equations (j|) and (||) were obtained starting from N coefficients of the power series S[t). They are 
therefore satisfied a priori only modulo . We compute more terms of S{t), say 2N, and check whether the 
same equations still hold modulo t^^ . If this is the case, chances increase that the guessed equations also hold 
at infinite precision. 

2.4.3 Analytic sieve: Singularity analysis 

By Theorem^ the minimal order operators for power series hke S{t) — F{t; 0, 0) and S{t) — F{t; 1,1) must 
have only regular singularities (including the point at infinity) and their exponents must be rational numbers. 

Example 6 (continued) The differential operator Ci^i is Fuchsian. Indeed, a (fully automated) local sin- 
gularity analysis shows that the set of its singular points { — 1, 0, oo, ^, |, —^(1 i iVS)} is formed solely 
of regular singularities. Moreover, the indicial polynomials of Ci^i are, respectively: t{t — l){t — 2)(2t — 
1), t(t - l)(2t + l){t + 1), [t -5){t- l){t -2){t-4),{t + 1)1(41 - l)(4t + 1), t(t -l)(t- 2)(t - 4), and 
t{t — 2){2t — 3)(t — 1). Their roots are the rational exponents of the singularities. 
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2.4.4 Arithmetic sieve: G -series and global nilpotence 

Last, but not least, we check an arithmetic property of the guessed differential equations by exploiting the fact 
that those expected to arise in our combinatorial context are very special. 

Indeed, by a theorem due to the Chudnovsky brothers the minimal order differential operator £ € 
Q[t]{Dt) killing a G-series enjoys a remarkable arithmetic property: C is globally nilpotent. By definition, 
this means that for almost every prime number p (i.e., for all with finitely many exceptions), there exists an 
integer /i > 1 such that the remainder of the Euclidean (right) division of D^^ by £ is congruent to zero 
modulo p [|l],|l6|]. 

From a computational view-point, a fine feature is that the nilpotence modulo p is checkable. If r denotes 
the order of C, let Mp be the p-curvature matrix of C, defined as the r x r matrix with entries in Q{t) whose 
entry is the coefficient of -Dj ^ in the remainder of the Euclidean (right) division of D^^^^^ by C Then, 
£ is nilpotent modulo p if and only if the matrix Mp is nilpotent modulo p [|l6[ 

In combination with Theorem |[ this yields a fast algorithmic filter: as soon as we guess a candidate 
differential equation satisfied by a generating series which is suspected to be a G-series (e.g. by 1, 1)), 
we check whether its p-curvature is nilpotent, say modulo the first 50 primes for which the reduced operator 
£ mod p is well-defined. If the p-curvature matrix of £ is nilpotent modulo p for all those primes p, then the 
guessed equation is, with very high probability, the correct one. 

We push even further this arithmetic sieving. A famous conjecture, attributed to Grothendieck, asserts that 
the differential equation £{S) — possesses a basis of algebraic solutions (over Q{x)) if and only if its p- 
curvature matrix Mp is zero modulo p for almost all primes p. Even if the conjecture is, for the moment, fully 
proved only for order one operators and partially in the other cases [|r3|, we freely use it as an oracle to detect 
whether a guessed differential equation has a basis of algebraic solutions. For instance, the computation of the 
p-curvature of an order 1 1 differential operator with polynomial coefficients of degree 96 in t, was one of the 
key points in our discovery [|| that the trivariate generating function for Gessel walks is algebraic. 

Example 7 (continued) The 5-curvature matrix Mz{t) of the differential operator £11 in (^ has the form 
-^M^{t), where d{t) = [it - iyt^{t+ \ f{9t^ + 3t + l)'5(3i - 4) and Mz{t) is a A x A matrix with 
polynomial entries in Q[/;] of degree at most 27. The characteristic polynomial XM^ of M^ reads 

+ ^ m{t) {3t - ir + ^ m) (3i - 1)'^ + N,{t) T + Noitl 

where Nq, Ni, N2, N3 are irreducible polynomials in Z[f], of degree, respectively, 21,26,26,21 and with 
coefficients having at most 20 decimal digits. 

The polynomial Xi^h obviously equals modulo p ~ 5, so the 5-curvature of £1,1 is nilpotent (but 
not zerop^ modulo 5. In fact, for all the primes 7 < p < 100, the p-curvature matrix of £1,1 is also 
nilpotent modulo p; it is even zero modulo p. Under the assumption that Grothendieck' s conjecture is true, 
this indicates that £11 admits a basis of algebraic solutions, and so provides independent evidence that also 
S{t) = F{t; 1, 1) is algebraic. 

3 Empirical Results in 2D 

In this section, we consider the total number of walks only, i.e., the generating function F{t] 1, 1). Because 
of symmetries, the 256 possible step sets give rise to 92 different sequences only. By inspection of the first 
N — 1000 terms, we found that 36 of them appear to be D-finite: 19 are algebraic and 17 are transcendental. 
The D-finite step sets, together with the sizes of the equations we discovered, are listed in Table [l| in the 
appendix. (There, and below, step sets are represented by compact pictograms, e.g. : for © = /*, [}.) 



(ill) Modulo 5, the curvature matrix M^(t) has as minimal polynomial. 
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3. 1 Combinatorial Observations 

Our classification matches the results of Bousquet-Melou and Mishna [^: for every sequence they prove D- 
finite our software found a recurrence and a differential equation, and whenever a series is algebraic indeed, our 
programs recognized it. Moreover, we found no recurrence or differential equation for any step set conjectured 
non-D-finite by Bousquet-Melou and Mishna. This strengthens the evidence in favor of the conjectured non- 
D-finiteness of these cases. 

3.2 Algebraic Observations 

All but two of the minimal polynomials of the algebraic series share the property that they define a curve of 
genus 0. As a consequence, there exists a rational parametrization in all these cases. For example, for the 
Kreweras step set : , the minimal polynomial Pi.i given in defines a curve parameterized by 

T(u) = \ and t{u) = , , , 

^ ' (u + 9){m2 + 24u+147) ^ ' a{u)' 

where a{u) = (u6+66u5+1827u''+27180w3+229431u2+1042866w+1995717)/(m+11)(u2+22u+125)2, 
i.e., for these rational functions we have 

Fi,i(r(u),t(w)) = 0. 

The two algebraic series that do not admit a rational parametrization belong to the step sets , . • (reverse 
Kreweras) and (Gessel's). Their genus is 1. 

Another feature of the series which we found to be algebraic is that they all admit closed forms in terms of 
(nested) radical expressions. For example, for the Kreweras step set, we find that F{t; 1, 1) is equal to 

i - V3) (216^3 + 1) (t - M^f - 2it{m^ - 15t + l)a(t) + {i + Vs) a{tY 

6it3(3i _ l)3a(i) 




where i = V^anda(t) = ^ 24:^J 3t^ {3t ~ 1)^ (91"^ + 3t + if ~ {3t - 1)3 (5832i6 + 540^^ - 1). Such 
representations can be found by appealing to the built-in equation solvers of Maple and Mathematica applied 
to the equation Pi i = 0. Both features are remarkable because, among all algebraic power series, only a few 
are rationally parameterizable or expressible in terms of radicals. 

Also the transcendental D-finite series appear to have some special properties. Being D-finite, these series 
are annihilated by some linear differential operator 

C = co{t) + Ci{t)Dt + ■■■ + Cr{t)Dl e Q[t](A). 

According to the DFactor command from Maple's DEtools package, all the operators can be factorized into a 
product of one irreducible operator of order 2 and several operators of order 1 . As all the operators are globally 
nilpotent, so are all their factors 1 16 1^ . 



We can therefore expect that every solution of these factors can be written as a sum of terms of the form 



R{ty ■ 2F1 



a 13 
7 



Z{t)\, (7) 



where R and Z are rational functions in Q(i) and q, /3, 7, 5 are rational numbers. Indeed, Dwork [|l6| Item 
7.4] has conjectured that any globally nilpotent second order differential equation has either algebraic solu- 
tions or is gauge equivalent to a weak puUback of a Gauss hypergeometric differential equation with rational 
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parameters. This conjecture was disproved by Krammer [ [Zq ] and recently by Dettweiler and Reiter 
the counter-examples given in these papers require involved tools in algebraic geometry (arithmetic triangle 
groups, systems associated to periods of Shimura curves, . . . ) 

We are therefore in a win-win situation: either the second order operators appearing as factors of our oper- 
ators admit only solutions which are indeed sums of terms of the form (^, or there is a simple combinatorial 
counter-example to Dwork's conjecture. Let us illustrate this on one of the most simple examples, the step set 
; . ; . We find here the differential operator 



4(32t^ - m - 1) + 4(8i - l){20t^ -M - 1)A + t(4t - l)(112i^ - 5)Dt + r (4i - iy{At + l)Di 

which Maple factors into 

(2(192t^ - 56t^ - 6i + 1) + 4(24^^ - l)(4t - l)tDt + {At - if {At + l)t^Dl) (l/t + A)- 

With the help of Maple's built-in differential equation solver (the dsolve command), it can be found that the 
differential operator gives rise to the representation 



(Incidentally, this solution can also be expressed in terms of elliptic functions.) We believe that all the tran- 
scendental D-finite generating functions for any step set admit a representation as (a nested integral of) such 
an expression. The solvers of Maple and Mathematica, however, are able to discover such a representation 
only in the simplest cases. (Note that at present, no complete algorithm is known that is capable of finding 
general pullback representations.) 

3.3 Analytic Observations 

By Theorem H, all the coefficient sequences grow like Kn°'p'"' log(rt)'' for some constants k, p, a, (3 (we only 
care about the dominant part of their asymptotic expansions). From the differential equation or the recurrence 
equation, we can determine p, a, and (3 exactly as roots of characteristic polynomials and indicial equations, 
respectively. (See [|3[ |l^] on how this is done.) We find that /3 = in all cases. Knowing the recurrence, 
we can also compute easily tens of thousands of sequence terms. With the help of convergence acceleration 
techniques [|l2| applied to so many terms, it is possible to determine the remaining constant n to an accuracy 
of thirty digits or more. With that many digits, it makes sense to search systematically for potential exact 
expressions of these constants using Plouffe's inverter [|o| and/or algorithms like LLL and PSLQ [^]. We 
actually found "closed form" expressions for all these constants. They are included in Table |l] in the appendix. 

By Theorem |lj the numbers p are bounded by the cardinality of the step set 6. It turns out that p = |S| 
unless the vector sum of the elements of the step set points outside the first quadrant. In these cases, p is an 
algebraic number of degree 2 (e.g., p = 1 + 2^/2 for the step set • . : ). For a, we found only non-positive 
numbers. Note that a being a negative integer implies that the corresponding series is transcendental |p^]. 

All the constants k have the form up'^°(t)^i (t>^ ■ ■ ■ 0^'', where the (j>i are usually small integers, the are 
rational numbers, and u is I/tt if F{t] 1, 1) is transcendental, and l/r(a + 1) if F{t\ 1, 1) is algebraic. 

There are some cases where the 0^ are not integers. Among them, very strange is only the case of the step 
set l\\ , for which we found r = 1, eo = 7/2, ei = 1/2, p = 2 + 2^6 and (j)i = (1137 + 468\/6) /152000. 
This last number may look like a guessing artefact at first glance, but we trust in its correctness, because the 
number of correct digits exceeds by far the number of correct digits to be expected from an artefact. 

4 Empirical Results in 3D 

We have investigated walks in three dimensions confined to the first octant with step sets of up to five elements. 
A priori, there are 83682 such step sets, and they give rise to 3334 different sequences. Of those, we have 
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computed the first N = 400 terms of the generating function G{t] 1, 1, 1) of general walks, and searched for 
potential differential equations, algebraic equations, and recurrence equations. We found that 134 sequences 
appear to be D-finite, and among those, 50 appear to be algebraic. 

4. 1 Combinatorial Observations 

For some of the sequences, it can be realized that their D-finiteness or algebraicity is a consequence of the 
D-finiteness or algebraicity of a certain 2D walk. For example, the sequence corresponding to the step set 

::: 1,4,17,75,339,1558,7247,34016,160795,764388,... (A026378) 

is readily seen to be D-finite, since it may be regarded as a variation of the 2D step set • . • in which the 
step t appears in two copies and empty steps are allowed. (Here and below, a three dimensional step set is 
depicted in three separate slices: first the arrows tops of the forms (x, y, —1), then (a;, y, 0), then (x, y, 1). For 
example, the step set above is {(—1, 0, 0), (0,1,0), (1,0,0), (0,0,1), (0, 1, 1)}. The given numbers are the 
first coefficients in the expansion of G{t\ 1, 1, 1).) 

Discarding these cases from consideration, we are left with 35 different sequences whose generating series 
appear to be D-finite; among those, three appear to be algebraic. Their step sets are given in the appendix. 

We were not able to find an equation for the step set 

\ \\ 1,1,4,7,28,70,280,787,3148,9526,38104,... (A149080) 

which is symmetric about all three axes, not even with 800 terms instead of 400. Also the step set 

\\' 1,1,4,13,40,136,496,1753,6256,22912,85216,... (A149424) 

which enjoys a rotational symmetry about the middle line of the first octant, and which may be viewed as a 
three dimensional analogue of Kreweras's step set, appears to be non-D-finite, even when 800 terms are taken 
into account. 

For walks in the quarter plane, it is conjectured in Section 3] that D-finiteness is preserved under re- 
versing arrows, i.e., the generating function for a step set 6 is D-finite if and only if the generating function for 
the step set &' is, when &' is obtained from S by reversing all arrows. Our computations do not suggest that 
this criterion also applies in 3D. Among the 134 sequences we found D-finite, there are 42 which correspond 
to step sets in 6 for whose counterpart in & we were not able to find an equation. Among those, there are 
some which satisfy only very large equations, so that chances are that they remain D-finite upon reversing 
arrows, but with equations which are too large for us to find. Others satisfy quite small equations, for example 
the sequence A026378 whose step set is given above. 

4.2 Algebraic Observations 

As in the 2D case, it turns out that most of the minimal polynomials of the algebraic series define curves of 
genus 0, which therefore can be rationally parameterized. There are twelve cases of genus 1, these are elliptic 
curves. Some of them turn out to be isomorphic (over Q). For example, those corresponding to the cases 

::: 1,1,4,11,32,110,360,1163,4112,14066,47848,... (A149232) 

::: •/ u\ 1,2,7,27,105,426,1787,7590,32633,142152,624659,... (A150591) 

::: ./ v' 1,2,10,40,176,808,3720,17152,81440,384448,... (A151023) 
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all have 1728 as j-invariant. They originate from the 2D Kreweras walks. Most interestingly, there are also 
three step sets originating from the 2D reverse Kreweras walks ( , . • ) for which the genus is 5 (!). 

For the transcendental series, we could observe the same phenomenon as in 2D: all the operators factor as 
a product of a single irreducible operator of order two and several operators of order one. We therefore expect 
again that all these series admit a representation as a hypergeometric pullback. As an example, the generating 
function G(t; 1, 1, 1) of the sequence 

-I 1,1,2,4,10,25,70,196,588,1764,... (A005817) 
can be written in the form 

4*2-2^+1 9, 9 , 9., -i 9 s A/2 1/2 



48i3 



■ x-^{4:X^ -2x + iy^(^{l - 16x^){'i8x^ - 12x^ - 8x - l)2i^i (^^"^ ^ 



16x2 



3 ..2 ^^-1/2 1/2 



(112x^-4x^-8x- 1)2^1 



1 



16x^ ]dx. 



This representation was found by Mark van Hoeij. It is beyond the scope of the standard tools of Maple or 
Mathematica. 

4.3 Analytic Observations 

Also concerning asymptotics, similar remarks apply as in 2D. All coefficient sequences grow like Kn°'p" for 
some constants k, a, p, where p is an integer or an algebraic number of degree 2 and a is a non-positive 
number We have not gone through the laborious task of determining the constants k. 

Acknowledgments. We are grateful to the anonymous referees for their pertinent suggestions which helped 
to improve the presentation of this paper. We thank Pierre Nicodeme and Bruno Salvy, who carefully read 
a preliminary version and made several useful comments. We also thank Mark van Hoeij for computing the 



example pullback in 3D given in Section 4.2. 
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Appendix 

Table 1 D-finite series and their step sets in 2D. The equation sizes columns refer to (minimal) recurrence 
equation, differential equation, and algebraic equation, respectively. Example: The series F{t; 1, 1) for Krew- 
eras walks (A151265) satisfies a differential equation of order 4 with polynomial coefficients of degree 9 and 
an algebraic equation P{F{t; 1, 1), i) = for a polynomial P{T, t) of degree 6 in T and 8 in t. The coefficient 
sequence of F{t; 1, 1) satisfies a recurrence equation of order 6 with polynomial coefficients of degree 4. The 
labels used in the columns "OEIS Tag" are taken from Sloane's On-Line Encyclopedia of Integer Sequences 
http : //www . research . att . com/~n jas/ sequences/ , Constants in the asymptotics columns are 

abbreviated A = l + V2,B = 1 + 2^2, C = l + \/3,L» = l + 2V3,E = 76(379 + 156^6) (!),F = 1 + ^6. 



Steps 



Equation sizes 



Asymptotics 



OEIS Tag 



Steps 



Equation sizes 



Asymptotics 



1, 

2, 1 

2, 1 

3, 1 
6,4 
7,4 

2,3 
2, 3 

2, 3 

3, 1 
3, 1 
9, 18 

3, 1 

8, 15 
2, 1 

4, 5 
7, 14 

9, 18 



1, 1 

2,3 

2,3 
2,5 
4,9 

4, 12 

3,5 
3,5 

3,5 
2,5 
2,5 

5, 24 
2,5 

5, 20 

2,3 

3, 8 

5, 18 

5, 24 



1, 1 

2,2 

2,2 
2,2 
6,8 
6,8 



8,9 
2,2 
2,2 

2,2 
2,2 



V2 2" 

r(I)^ 

3^3 3" 

2r(i)^ 

12v^ 23"/2 
r(i) n3/2 
2V2 3" 

3^3 3" 

8 4" 
TT 

2 4" 
7r n 
4^3 4" 
3r(|) 
1 

_4" 
2 

^2^3/2 

23/4r(i)^ 

12Vm (^24)" 
4/3^.2^)3/2 j^n 



8r(|) 

V5 



l3/2 



2V2r(|) 
3v^ 6" 

2r(|) 

2^3 6" 

3r(|)7s 

^7/3 7" 

3r(|)7^ 



A00()079 
A00()244 
A005773 
A151255 
A151266 
AI51281 

A005566 
A06()899 

A128386 
A151261 
A151291 
A151287 
A151302 

A151318 
A151297 
A151323 
A151314 
A151331 





1 

8 

10 
, 1 



3, 1 

5, 8 

6, 10 

7, 11 
9, 18 

2, 1 

7, 11 

2, 1 

9, 18 

3,4 



1, 1 
1, 1 
2,3 

4, 16 

5, 16 
2,5 
3,4 
2,3 
2,5 

4, 15 

5, 15 

5, 19 
5,24 

2,3 

5, 18 

2,3 

5,24 

3,6 



1, 1 

1, 1 

2, 2 



2, 2 



2, 2 



2, 2 



2, 2 



4,4 



2" 



V3 3" 

24^2 23"/2 
7r ?i2 
V3 3" 

2r(i)7^ 

2 

4 4" 

TT n 
s/2 4" 

6v^ 2"3"''2 

12^3 2"3"/2 



3r(i) 
VSA"^/^ (2A)" 
7r n2 
V5 5" 

^5/2 5" 

V3Cy^ (2C)" 

27r n2 
V233/4 6" 

_BF^/2 (2F)" 
5V957r n2 
8 8" 
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Table 2 Conjecturally algebraic series and their step sets in 3D. Step set figures are as in Section ^ Equation 
sizes are as m Table |l[ 



First terms (OEIS Tag) 


Step sets 


Equation sizes 


1, 1,4, 10, 37, 121,451, 1639, .. 


(A025237) 


... 

:•: :.: ::: 


■•■ :• ::: 

... • . . 


2, 1 


2,3 


2,2 


1, 1,5, 15,51, 199, 755,2789, .. 


(A149576) 


... • . . • • • 
... • . ■ ■ ■ 
... • . . 


... ... ■ • • 

... • ■ ■ ■ 
... • . . 


11, 22 


7, 31 


12, 17 


1, 2, 4, 14, 46, 134, 502, 1820, . . 


(A149847) 


::: :■• 

• ■ ■ 


... . • 


8,6 


4, 16 


6,9 



Table 3 Conjecturally transcendental D-finite generating series and their step sets in 3D. The equation sizes 
columns refer to (minimal) recurrence equations, and differential equations, respectively. 



OEIS Tag 


Step sets 


Equation sizes 


OEIS Tag 


Step sets 


Equation sizes 


A 1 A9.C\f\C\ 


• - - ■ ■ • 

• ■ 

• ■ ■ • ■ 

• - ■ ■ ■ • 

• ■ ■ 

• ■ 

• ■ • ' ' ' 

• ■ 

• ■ ■ 


• • ■ 

• - - 

• • ■ 


9, 17 


5,28 






• 

• 
• 


• - - 
• 

• ■ ■ 

• - ■ 
• 

• ■ ■ 

• ■ ■ 
• 

• ■ ■ 


- • - 

- • - 

• • • 


7, 10 


5, 17 




• - • ' ' ' 

• ■ ■ 

• ■ ■ ■ ■ ■ 


• • - 






A 149090 


• 




- ■ • 


• ■ - 

• ■ - 

• ■ - 


9, 17 


5,28 


A 1495 89 


■ • ■ 


• • • 

• • ■ 

• • ■ 

• - • 

• • ■ 

• - - 


10,21 


6, 29 




• 
• 


• 
• 
• 


- ■ • 


• • • 

• ■ - 

• ■ • 

• ■ - 






A005817 


... ■ • ■ 


- • - 


2,2 


3,4 


A 148005 


• 




• - ■ 


• ■ - 


5, 8 


4, 15 


A148052 


• - • • ■ ■ 

... • - ■ 


• - - 


7, 18 


6, 27 


A148068 


• 
• 






• ■ - 

• ■ - 


7, 17 


6, 25 


A148072 


... • ■ ■ 
... • - ■ 




12,57 


10, 69 


A148162 








• ■ • 

• ■ • 


4,3 


3,6 


A148284 


... • • ■ 
■ ■ • * 

... • - ■ 




14,57 


10,71 


A148331 


• 


• 


• ■ • 




11,43 


9,53 


A148507 


• ■ ■ ■ ■ ■ 
... • • 

... • • ■ 


• - - 


4,6 


4, 11 


A 148525 






• ■ ■ 

• 

• ■ ■ 




7, 16 


6, 25 


A148548 


... • ■ ■ 
• - ■ • ■ ■ 


• - - 


7, 19 


6, 28 


A148689 




• 


• ■ ■ 

• ■ ■ 




8, 25 


8,31 


A148703 




• • ■ 

• - ■ 

• - ■ 


4,3 


3,6 


A 148790 




• 


• 




6, 12 


5, 18 


A148934 


■ • ■ 

... • - • 


• • • 


5,5 


4, 11 


A149279 




• 


• - ■ 
• 

• - ■ 




14, 62 


10, 75 


A149290 


... • - • 


• • • 


11,53 


9,61 


A149363 




• 


- - • 
• 

• ■ • 


• ■ • 


7, 16 


6, 24 


A149632 


■ ■ • • ■ ■ 

... • 
... • ■ ■ 


• • • 


7, 11 


5, 16 


A149713 








• ■ • 

• ■ • 


8, 22 


7,29 


A150054 


... • - ■ 
■ • - • 

... • • ■ 




12, 39 


9, 52 


A150370 






• ■ ■ 

• 

• ■ ■ 


- ■ • 


14, 62 


10, 75 


A150410 


■ ■ • 

... • • 

... - • 


• • • 


4,6 


4, 11 


A 150471 






- - • 

- • ■ 




12, 33 


8,42 


A 150499 


■ • ■ 

... • • 

... ■ • ■ 


- • - 


14, 48 


9,61 


A150764 






• ■ • 

• ■ • 


- ■ • 


7, 13 


6, 19 


A150950 


■ • ■ 


- - • 

• • • 

• • ■ 


8, 23 


7, 29 


A151053 






• • ■ 

• 

• - ■ 




14, 38 


9, 48 



